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In this work, a physical system described by Hamiltonian = Ho + Vtj(x,t) consisted of 
a solvable model Hq and external random and time-dependent potential V„ (x, t) is investigated. 

Under the conditions that the average external potential with respect to the configuration uj is 
constant in time, and, for each configuration, the potential changes smoothly that the evolution of 
the system follows Schrodinger dynamics, the mean-dynamics can be derived from taking average 
of the equation with respect to configuration parameter uj. It provides extra contributions from the 
deviations of the Hamiltonian and evolved state along the time to the Heisenberg and Liouville- 
von Neumann equations. Consequently, the Kubo’s formula and the fluctuation-dissipation relation 
obtained from the construction is modified in the sense that the contribution from the information 
of randomness and memory effect from time-dependence are present. 
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I. INTRODUCTION 


The fluctuation-dissipation relation (FDR) is one of the most well-known formulae in statistical physics especially 
in weakly perturbed systems around equilibria. The relation was first coined by H. B. Callen and T. A. Welton in 
1951 [ij and developed by R. Kubo in 1957 Q. It states that the rate of energy dissipation, which one can obtain 
from a measurement, can be described via a fluctuation of the system in terms of the correlation amongst the group 
of considered quantities It is applicable to explain many physical systems, especially in a measurement process, 

since one can perturb the considered system to obtain a response function and infer properties of the system. 

However, reports over the decades have indicated that the FDR does not hold appropriately in many situations 
such as glassy or driven systems. In a system with very slow change toward equilibrium, which can be due to the 
composition of many complicate effects concerning additional interactions, the FDR could be modified into Quasi- 
FDR that the correlation and response functions are extended into a more general form ii- Such extension covers, 
both theoretically and experimentally many cases of models which concern, for instances, the long-time relaxation, 
structural glasses [EHLspin glasses Ising model with dipolar interactions 0 , spin-boson model [l^ , and the 

Glauber-Ising chain [lfi| . Also, a driven system is another type of Hamiltonian model which is distracted from the 
formal FDR by a different reason: time dependence of the Hamiltonian. Nonetheless, such a system provides similar 
characters of dynamics as the glassy class. That is, the significant slowness of evolution toward equilibrium or toward 
steady state. In essence the time dependent Hamiltonian effect can arise from the memory content of the dynamics 
which prevents the system to lie in or near the equilibrium component [itI - I^ . 

Indeed, in literature, the derivations of such deviation of FDR and Kubo’s formula can be obtain by modifying 
the standard approach involved with perturbation from a bath called linear resMnse formalism Q. There are many 
suggestions of the modifications presented which the readers can consult Ref.jl, [23l - l^ and references therein for 
further information. Almost all of them gave similar revisions of the modified-FDR; namely, the standard one but 
added by an extra term into the equation. In this work, we introduce an alternative explanation to obtain a modified 
Kubo’s formula and the modified-FDR by using the role of random, time-dependent potential into consideration. 
Also, the discussions on characteristics of the model and related topics will be given. 


II. FORMULATION 

A. A Model 

We consider a system described by a time-dependent Hamiltonian 

H<.(x,t) = Ho + V^(x,t) (1) 

(denoted by Htj(<) or for short), where Hg is a fixed Hamiltonian, and V^(x, t) is multiplicative operator by 
a random smooth function ?;^(x, t) of t; that is V^(x, t)(^ = ?;^(x, t)i^ for each ip in the Hilbert space %. Here the 
parameter w denotes a configuration in a probability space (D,/i), with a measure /i is a measure associated with 
the configurations oj. It can be assumed that Hg is solvable so that one can formally describe all properties of the 
system ^g. V(^(x, t) describes effects from an external system Suj connecting to the system S'g, but process random 
fluctuation. The role of Suj is to supply or extract energy by random amount of at time t. This causes the 

measured total energy of the composite system, denoted by 5' = ^g U 5^^, to depend on time t. In our consideration, 
the composite system of the previous systems is also weakly contacted with a reservoir B with a bath parameter 
j3 = l/ksT, where T is a bath temperature to make it more realistic in thermodynamic sense, and more convenient 
to calculate statistical physical quantities. Remark that a glassy-type system or a driven system can be described in 
this framework by assigning a specific random profile or time-dependence profile to the potential (t). 

Suppose {(jik) is a set of eigenfunctions of Hg in the Hilbert space % corresponding to eigenvalues {A^}, where 
k = l,2,...,dg := dim'H. We will restrict our consideration to only the Hilbert space TL. However, one can also 
explain the composite system S by another Hilbert space which is possibly larger or smaller than 'H because of the 
additional term of Hamiltonian V^(x, t). Note that the interaction with the bath is sufficiently weak that it does not 
change the structure of the composite system, or, in our case, the system Sq (and H). For the latter, we consider the 
case that the Hamiltonian Hg is that of a dg-level system. For such a system, one can obtain an equilibrium state or 
Gibbs state which maximizes the von Neumann entropy among all states with a given energy Eq [23, [23 as 

■- Tr[e-/3Ho] - ^ J 


\4’k){4’k\ 


( 2 ) 
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which is equilibrium for only the case Ui^(x, t) = 0 at all time t. In this case, Eq = := Tr[cr^H], where 

(A)p = Tr(pA) is an average over a certain state p of the operator A. On the contrary, for the system governed 
by H,,, it becomes more complicate to find an equilibrium state from a given energy because the Hamiltonian 
is time-dependent; hence the energy Tr (cr^Htj(<)) of the system in the state ap is no longer a constant of motion. 

In this scenario, for a path w, the state of the form - 3 = 5 — is not an equilibrium state. Moreover, since also 

depends on configuration w, one may not obtain the same result from the different preparations or measurements. 
Therefore, we will impose additional assumptions to define an appropriate constant of motion to develop another 
version of equilibrium state for overall evolution including the effect of anomalous potential V^(x, t). 


B. Equilibrium State 

In order to construct an equilibrium state, we add additional assumptions to the external system Sui- For each 
configuration w G H, Vui{x,t) is supposed to be a smooth function of t. We denote by (Xi^) := J^X^dp.{uj) the 
expectation or w-average of random variable X^, with respect to the measure p, on the probability space (H,/i). For 
simplicity, we assume that: 

V := (V^(x, t)) is a constant (multiplicative operator) in time. (3) 

From the condition (jS)), we find that Hq and V are time-independent and are also the constants of motion. Thus, an 
equilibrium state in this case, called adjusted equilibrium state, can be written as 

g-/3(Ho-Si^V) 

= Xi-(e-/5(Ho-ti^V))’ 

where /3 and v are Lagrange multipliers arisen from the formulation in Ref.jl^. Without loss of generality, we set 
= I by absorbing the scaling into the random part of Hamiltonian. Note that the parameter /3 must equal with the 
inverse temperature of the bath when V^(x, t) = 0 for all time t. Moreover, by employing the equivalence between 
the Schrodinger and Heisenberg pictures, for an evolution of the operator A A(t) and its dual evolution for a state 
p p{t), one obtains 


(A)p(t) = Tr(p(<)A) = Tr(pA(t)) = (A(t))p. 


Exact forms of the evolution maps will be discussed in the following subsection. 


( 5 ) 


C. Dynamics of Operators and States 

For a fixed configuration w, the dynamics of the system in this specific realization is well defined formally since 
Uij(x, t) is assumed to be a smooth function. The evolution of a self adjoint operator observable A and its dual 
dynamics follow 

j^A^it) = iCl{A^{t)) :=z[H^,A„(t)]; ( 6 ) 

= -iCl{p^{t)) = [H^,p<^(t)] . (7) 

Furthermore, one can also define the dynamics generated by a mean Hamiltonian H = Hq -I- V. The mean dynamics 
of operator A and p can be obtained by taking average over all configurations in Eq.® and Eq.([7]). Thus, 

^A{t)=i[n,A{t)]+iCt (H,A)^ 
at 

^p{t) = -i [H, p{t)] - iCt (H, p ), 

where A{t) and p{t) denote the w-average of A^(<) and Pui{t), respectively; and 

Ct(H,-) :=E,,{[SllUt),UUt)i-)]) (10) 

with (5H(^(t) := H^(t) —H and ^^^(t) being unitary evolution map from the initial time t = 0 to time t corresponding to 
Eq.© for a specific configuration w, is a non-homogeneous contribution from the mean-deviation of the configurations 


( 8 ) 

( 9 ) 
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about the mean Hamiltonian H. One can see that the argument in Eq. m is a Schrodinger operator and also an 
initial state of the evolution in the previous system of equations Eqs. dH)-©. Moreover, from Eq.©, one can see 
that the term C* (H, p) = Ct (V, p) contains the information of dependence between the Hamiltonian and the density 
operator or the state. Ct (H, •) vanishes when commute with and Pu;{t) at any time t. 


III. MODIFIED KUDO’S FORMULA AND FLUCTUATION AND DISSIPATION THEOREM 

A. Effects from the Bath 

From Eq. one obtains an ensemble average of a physical quantity corresponding to the operator A by employing 
the equivalence between the Heisenberg and Schrodinger pictures, so one need investigate only the dynamics of the 
state p. Here we will consider the linear response formulation as define in Ref.[^ to find the behaviour of the system 
when it is weakly contacted with a bath. First of all, we define a small interaction term with the bath as 


4 (-) := [B,-] 


( 11 ) 


where B is an operator representing the effects from the bath, and /ib is its corresponding c-number. Adding this 
term to Eq.® leads to a perturbed evolution[2^ as 


—p(t) = -i{L- L\) p{t) - iCt (H, p ), 


( 12 ) 


where C is the Liouville operator corresponding to the mean Hamiltonian H. Its formal solution can be written as 


p(t) = e-**^p(0) + t f 

+ r]{t) + i [ Cjr]{s)ds 

Jo 


+0{h 


Bh 


where 


Jo 


(13) 


(14) 


In the case, the initial state p{0) is an adjusted equilibrium state a'^, so it follows that e = cr'^. Combining 

Eas. lfT^ and ® and performing iteration to calculate the average of the observable A; 


(A(t))p = (A)^ + i f (Tr {a'p [A(t), B(s)]) ^ 3 ( 5 )) ds 
Jo 

p{t)A + i [ e-*(*-")^£f77(s)Ads 

Jo 


+Tr 


+Oihl 


(15) 


From the expression above, one can see that an average of physical quantity is deviated not only by a perturbation 
from a bath solely (the second term on the right hand side), but also by the non-homogeneous term r]{t) in the third 
term and coupling between the bath and random effects in the forth term. In the sense that the bath is a measurement 
probe, a linear response which is related to a measured quantity will be deviated by those effects also. 


B. Modified Kubo’s formula and Fluctuation Dissipation Relation 

The linear response function is defined by 




d{A{t)), 
dh-B {t') 


h ^—0 


( 16 ) 
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Then we obtain, for 0 < t' < t, 

XAB{t, t') = 2i9{t - {a'p [A{t),B{t')]) + Aab(<, t')}, 

where 

pt' 


(17) 


1 , 


AAB{t,t') := -Tr 




(H,p) 


ds] A 


(18) 


In physical sense, the term Tr [A(t), B(t')]^ can be related to an symmetric correlation Q, where correlation is 
usually defined as a two-epoch time correlation function between operators A and B in a given state. In this case: 


Cab := {A{t)B{t'))p. 


(19) 


One can then define symmetric and anti-symmetric correlations respectively by 

CiB(t,t'):=([A(t),B(t')])/3, (20) 

Cl^{t,t')-.= {{A{t),B{t')})p. (21) 

where {A, B} := AB -f BA. Thus, the modified form of the Kubo’s formula can be written as 

XAB (i, t') = 2te{t - t') [C^B (i, t') + Aab (i, i')] ■ (22) 

Note that, from a typical set-up of linear response, /ib( 0 to b® zero for t < t' and then turned on at t' to reserve 
the causality of the measurement - the effect from the probe cannot influence the system before they are in contact 
at the time t = t' [laiiii. In other words, from the equation it can be said that the system is in adjusted equilibrium 
state at the initial time t = 0 and evolves until t = t' . Then, it makes weakly contact with the bath (which can well 
be a measurement probe), and they evolve together until the present time t. All information of the evolution of the 
system and of the concerned operator is therefore contained in Eq. (1221) . 

Now consider the conditions to obtain the original FDR from Ea. d^^ . Using the short-hand writing -A (t) 

when t' = 0, so that Cab {t, t') -A Cab (1), we obtain 

Cab {t) = Tr {a'pA{t)B) 

Tr(U(t)BUt(i/3)Ut(t)A) 

“ TrU(i^) 

Cab {t) = Cba {—t — 0) ■ (23) 

One can see that the relation above is another version of the Kubo-Martin-Schwinger(KMS) condition, and the 
evolution along imaginary time arises here. 

Suppose the interaction of the composite system with the bath begins at t' = 0 (or equivalently, the probe start at 
time t = 0 when the system is still in adjusted equilibrium state.) Then one can compute that Aab( 0 = 0 [s^; and 
Ea. (l22ll becomes 


XAB{t) = 2ie{t)C^^it). (24) 

Taking the Fourier transform of the last equation, one obtains 

Xab(A) = Cab (A) = (1 - Cab (A). (25) 

The last relation is a direct consequence of Eq. (1^^ . Furthermore, one finds that 

XAB (A) = tanh CXb (A) (26) 

which is similar to the standard fluctuation-dissipation relation in the sense that it is calculated from a Gibbs state 
but with a modification that the adjusted equilibrium state a'p contains the information of randomness. 
Nonetheless, notice that only the average (wrt w) of the random potential effects to the adjusted equilibrium state, 
resulting the shift to the measured energy. For the case that V = 0, despite having external w-fluctuation, the 
adjusted equilibrium become coincide with ap and the last relation is exactly the standard FDR. 
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IV. DISCUSSION AND CONCLUSION 

a. Connection to Other Works We revisit the model under consideration: (x, t) — Hq + (x, t), where Hg 

is a constant or fixed Hamiltonian, and Vtj(x, t) is multiplicative operator by a random smooth function t). The 
presence of a random external random potential introduces unknown degrees of freedom. Unlike the bath or reservoir 
model, all dynamical properties of the external system Si^ such as the relaxation time, dimension or volume are not 
specified - so that the system can evolve to equilibrium as contrary to the conventional formulation of bath. This 
idea is analogous to having unobserved degrees of freedom. Similar model was analysed by Budini et al in Ref. 

to derive the master equation of non-Markovian type. It suggests that the presence of additional degrees of freedom, 
random disturbance to the Hamiltonian of the composite system, provides the memory effect in the master equation 
( 3 ^ . In this work, instead of using the Hermitian operators to capture the unobserved degrees of freedom associated 
with the system with random coupling constant, the potential is employed in form of time-dependent random operator. 
One can formulate similar same scheme by setting that Vi^(x, t) = Ai^Q, where Q is a Hermitian operator, and Ai^ 
is a random real number representing the strength of the coupling, and by using the non-Markovian master equation 
instead of Eq.®. It can be seen that the potential can be reduced to time-independent, and the adjusted equilibrium 
state can be defined because the a w-average of the Hamiltonian yields a constant of motion in this case. The extension 
in the direction of the present work is still open in general, especially for physical interpretations and implementations. 

b. Connection to Stochastic Path Integral Formalism It is observed here that, as the commutator obeys the Lie 

structure, the term Ct (H,p) can be analysed as a derivative. In this instance, the deviated Hamiltonian := 

appearing in Ct (H, p) will generate a dynamics for an appropriate dynamical parameter, and the derivative 
(concerning a commutation) represents a change of quantities related to randomness. For example, let consider a naive 
form, for a particular configuration w, 


[6U^{t),puj{t) - p{t)]. 

Since a trace class Bi{H) = {A G B{'H) : Tr(A) < 00 }, where B{H) is an algebra of all bounded operator on "H, 
together with Hilbert-Schmidt inner product (A,B) = Tr [AIB] forms a separable Hilbert space H [13, P-33], clearly, 
the set of density operators is contained in H(H) or H, and any density operator p can be viewed as a vector in H 
while an action TZu:{-) := i [(5H^(t), •] is an operator on H. Let B{'H) denote a Banach space of all bounded operators 
on ?i, one can find that the inner product 

{p,T^u:ip)) = i (Trp[5H^(t),p]) = 0 

following the cyclic invariance of trace function. Thus, by Lumer-Phillips Theorem, the operator TZuji') is dissipative 
and can be a generator of contraction semigroup on Ft [Isl . Theorem 2.5]. However, the semigroup given here is 
determined up to the configuration w, so does a dynamical parameter (denoted by t^j) of the evolution. Therefore, 
we will obtain 

|^=7^.(p), (27) 

OTqj 

another governing equation for density operator p. Hence, Eq.® and Eq.(|n|), for the case of countable number of 
configurations case, become multi-dimensional problem (in parameters). When the cardinality of H is uncountable, 
the description of the problem in this direction will be considered as a version of path-integral or stochastic dynamics. 
Indeed, derivations in those fields can be adopted to analyse the problem in this scheme, and still remain open. 

c. Remark on Equilibrium State By the definition of the external potential alone, it can be seen that the dynamics 
of the system looks suspicious because it may not reach the equilibrium within a finite time. In fact, the equilibrium 
state, which maximizes the von Neumann entropy, is not well defined in this case because of lacking of information 
all over the time domain. Even though the Gibbs’ state can be defined for specific Hamiltonian at specific time, the 
time-dependence of the Hamiltonian causes ambiguity to the model by the fact that the Gibbs’ state defined here is 
dependent on time following the Hamiltonian. That says, when the Hamiltonian changes, the preferred direction of 
the dynamics in the state space also changes. Thus, the equilibrium condition depends on time and the dynamics as 
equilibration becomes questionable. In order to answer the question, one idea is that one can define an instantaneous 
equilibrium state which is given by the Hamiltonian at instant time and for specific configuration u. Then, the 
dynamics can be considered from a trajectory of such instantaneous equilibrium states, and the deviation from that 
state involves the entropy production along the path dynamics . On the contrary, in this work, by employing the 
assumption © that the energy corresponding to the w-average of Hamiltonian is a constant of motion, and one can 
define a single adjusted equilibrium state, irregarding to time, as in this formulation. Hence, the mean dynamics of 
the state can be derived from the equilibrium condition. 
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d. Remark on Entropy Production There remains interesting topic involving the entropy and information in views 
of the entropy function S{p) := —Tip In p. It is already known that the evolution of the system for a fixed configuration 
w is unitary, but the unknown random potential still affects the dynamics. Thus, one cannot know exactly in which 
the state the system lies. When we reduce the dynamics to the mean dynamics as in Eq.®, we expect entropy change 
by the averaging evolution process. To demonstrate this issue, consider Eq. 0, one can see that their entropies are 
unchanged along the time by invariance under the unitary transformation of the entropy function [^ . Because the 
entropy function is concave, it follows that 


S{p{t)) = S{E^ {p^t))) > E,, (Sip^t))) 

= E,, (5(p(0))) = 5(p(0)). 

Then, the margin of the entropies S{p{t)) — S{p{0)), which can be treated as the overall entropy production along the 
dynamics until the time t, can be positive. Additionally, for the case that Ct (H,/?) = 0 which is occur, for instant, 
when the random effect is absent; one can define an entropy production as 

Ent-B {t) ■- -iphs (i)Trp(t) [H, B] , (28) 

since the effect from the bath is treated as a perturbation and by choosing the adjusted equilibrium state as a 
reference state. Because the entropy production above is in linear order of /ib, the evolution is still close to the 
adiabatic regime. On the contrary, for the case that Ct (H,p) ^ 0, the evolution map corresponding to Eq.® may 
not admit a group or even a semi-group properties in general. The definition Ea. (|28l) is therefore not well defined for 
this case, and one can still expect a non-adiabatic effect in some situation. 

e. Remark on the Extra Term According to the derivation of the modified FDR in Section FlII B[ it is an advantage 
of using the mean dynamics Eq.®. As a result, it coincides with many previous works of the modified FDR: the 
correlation in the former FDR is replaced by a composition of the old one and another extra term. It was often 
suggested in literature that the function of the extra term is to explain the behaviour of the considered physical 
system out of the equilibrium regime. One idea is that the extra contribution to the response function is due to the 
dynamical activity, the time-symmetric part of action from the bath (or external system), which usually vanishes in 
equilibrium by the causality argument - the response must not occur before a measurement [s^. However, when the 
system is out of equilibrium or there are other parameters not included in the consideration, the dynamical activity 
seems to cause an amount of energy dissipation [23 - [^ . It has been also interpreted as total entropy production 

from the d yna mics, while the former term is treated as the entropy production of the medium or the considered 
system In our case, the extra term arises from the presence of 


Ct (H,p)=E^([m„(t),p,,(t)]) 


in the governing equation of the mean dynamics, where (5H,j(t) := H„(t) — H. One can see that it keeps all the 
information of the external random potential at any time through which the system evolves. Recall Eg. dTSl) for the 
explicit form of the extra term in the Kubo’s formula 








X 


(H,p) 


ds^ A 


(29) 


It can be seen that the actions along the dynamics are all included in the expression at the point they act. There 
are two evolutions therein: one maps from the initial time to the measurement time accumulating the effects from 
the non-homogeneous term Cg (H, p) interacted with a bath via the operator B along the time, and the second one 
)c system and bath evolve until the final time t. Therefore, it reflects the memory effect of the 

whole evolution, equipped with the non-homogeneous term from the random potential, in the measured quantity. 
In particular, the added random term disturbs the measured quantity to be deviate from a considered value in any 
experiment it involved with unless in the limit t' -A 0; making the state in the adjusted equilibrium state at the 
measurement time as shown in Ea. dMl) . 

/. Conclusion In summary, we find that the deviation of the FDR can be designed as an effect from a randomness. 
By considering a class of Schrodinger type equations indexed by configurations in a probability space and taking an 
configuration-expectation, the mean dynamics with non-homogeneous term will be obtained. A perturbation from a 
bath in linear order yield a deviation of average of a physical quantity which affects a linear response. Consequently, 
the modified-FDR is derived and, we find that it also relax to the standard FDR in the case that the probe begins to 
contact with the system when the latter has not yet reached an equilibrium, but is in an adjusted equilibrium. 
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Appendix A: Detailed Derivation 

For a fixed configuration w, the dynamics of the system in this specific realization w is well defined formally since 
Uij(x, t) is assumed to be a smooth function. Thus, the evolution can take the form 
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where T is a time-ordering operator with its reverse time-ordering operator denoted by T [^ : here, we set 0) := 

Uij (t) for simplicity. The evolution of the operator A follows 


with its dual dynamics 


■^A^{t) = iCl,{A^{t)) := i [H^^, A<^(t)], 

(Al) 

A^it) = So dsciA = ui(t)AU<.(t), 

(A2) 

= -iCl{p^{t)) := -i \R^,p^{t)] , 

(A3) 

Pu{t) = e~^'^So p = lJa,{t)pUl{t). 

(A4) 


However, we investigate its mean dynamics rather than the class of equations above, where the mean evolutions in 
Eqs. (|A1|) and (IA3I) are given by 

^A(^)=^[H,A(^)]+^Ct(H,A)^ (A5) 

^p{t) = -i [H, p{t)] - iCt (H, p) , (A6) 

where A{t) and p{t) denote the w-average of Ai^{t) and Pui{t), respectively; and 

Ct(H,-) :=E^{[SUUt)Mt){-)]) (A7) 

with := H,^(<)-H and U^{t){-) = being a unitary evolution map from the initial time t = 0 to 

time t corresponding to Ea. (IA3l) for a specific configuration uj. Indeed, Ct (H, •) can be treated as a non-homogeneous 
contribution from the mean-deviation of the configurations about the mean Hamiltonian H. For the Eo. ljABp which 
we will mainly focus has a formal solution as 

pit) = e"*‘^p(0) -b r]p{t), (A8) 

where 


Pp{t) := i f (H, p) ds. (A9) 

Jo 

One can see that, from the integral form above, the evolution map, in general, may not admit the semi-group properties 
by the presence of the non-homogeneous term Cg (H, p). 

Next, consider a perturbation from a bath C*j- = h^{t) [B, •]. The perturbed evolution can be written as 


dt 


p{t) = -i (£-4)p(t)-zCaH,p), 


(AlO) 


where C is the Liouville operator corresponding to the mean Hamiltonian H. Its formal solution can be written as 


p(t) = e-**^p(0) -b ^ 

Jo 


+r]pit) +i e' 
Jo 

+0{h^). 


-i{t-s)Cps 


C]fjp{s)ds 


(All) 


From the expression above, one can calculate the average of the observable A by applying the Schrddinger- 
Heisenberg equivalence and select p(0) = cr^, one obtains 

rt 


{A{t))p = {A)p + *Tr 


-bTr 




+0{hl). 

Next consider the integrand in the second term on right hand side of the equation, one can see that 


(A12) 
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Tr(e-*(*-^)^[B,4] a) =Tr(U(t)Ut(s)[B,a^]U(s)Ut(t)A) =Tr(a^[A(t),B(s)]), (A13) 

where U(t) = A{t) = U^(i)AU(i) and B(i) = U^(<)BU(t). The relations above arise from the cyclic invariance 

of trace operator and the fact that 

The linear response function is defined by 


Xab 


d{A{t)), 
dh-B {t') 


/iB—0 


Then we obtain, for 0 < t' < t, 


XATi{t,t') = 2ie{t - t'){Ti {a'l^ [A(t),B(t')]) + AAB(t,t')}) 

where 


(A14) 


(A15) 


AAB(t,t') := 



(H,p) 

ds] A 

\Jo 


) \ 


(A16) 


In physical sense, the term Tr [A(t), j can be related to an symmetric correlation [^; where correlation is 
usually defined as a two-epoch time correlation function between operators A and B in a given state. In this case: 


Cab := (A(t)B(t'))/3. 


(A17) 


One can then define symmetric and anti-symmetric correlations respectively by 

:=([A(t),B(t')])/3, (A18) 

C^bM :=({A(t),B(0})/3, (A19) 

where {A, B} := AB -|- BA. Thus, the modified form of the Kubo’s formula can be written as 

XAB [t, t') = 2i9{t - t') [CXb {t, t') + Aab (t, t')] , (A20) 

where 9{t) is a Heaviside function defined by d{t) = 1 for t > 0 and 6{t) = 0 for t < 0. Now consider the conditions 
to obtain the original FDR from Ea. (IA20(l . First, it follows that 

Cab (t,t') = Tr(4A(t)B(t')) 

= Tr {a'^A(t - t')B) 

= Cab (t - t', 0) (A21) 

by the cyclic invariance of trace. Thus, 

= (A22) 

Using the short-hand writing -A (t) when t' = 0, so that Cab -A Cab (t), we obtain 

Cab (t) = Tr {a'^A{t)B) 

Tr(U(t)BUt(i/3)Ut(t)A) 

“ TrU(i/3) 

= Tr (a^U'f(-t)Ul'(-t/3)BU(-i/3)U(-t)A) 

Cab (t) = Cba {—t — i/3) ■ (A23) 

One can see that the relation above is another version of the Kubo-Martin-Schwinger(KMS) condition, and the 
evolution along imaginary time arises here. 

Suppose the interaction of the composite system with the bath begins at t' = 0 (or equivalently, the probe start 
at time t' = 0 when the system is still in adjusted equilibrium state.) Then one can verify that AAB(t) = 0; and 
Ea. (IA20() becomes 


XAB(t) = 2ie{t)C^^{t). 


(A24) 
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It would be mentioned here that the resulting equation can be approached by another reason or assumption, for 
instance, when one forces the system so that Pui(t) always lies on its mean state p{t) - which is an alternative 
expression for mean-field approximation also. 

For an integrable function g{t), we define its Fourier transform as 

g{X) = T {g{t)) := -j= / dt, 

where A is a complex number. Taking the transform of Ea. (IA24l) . so 

Xab(A)=C'Xb(A) 

= (l_e-/5^)C'AB(A). (A25) 

The last relation is a direct consequence of Eq. (IA23|) . Furthermore, one finds that 

^AB(A) = (l + e-''^)CAB(A). (A26) 

Combine terms together in the equations above, one can obtain the standard fluctuation-dissipation relation 

Xab(A) = tanh (A27) 







